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Abstract 

Many-particle electron states in semiconductor quantum dots with carrier-mediated ferromag- 
netism are studied theoretically within the self-consistent Boltzmann equation formalism. De- 
pending on the conditions, a quantum dot may contain there phases: partially spin-polarized 
ferromagnetic, fully spin-polarized ferromagnetic, and paramagnetic phases. The physical prop- 
erties of many-body ferromagnetic confined systems come from the competing carrier-mediated 
ferromagnetic and Coulomb interactions. The magnetic phases in gated quantum dots with holes 
can be controlled by the voltage or via optical methods. 

PACS numbers: 73.21.La, 73.21.-b, 73.63.Kv, 78.67.Hc 
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I. INTRODUCTION 



Diluted magnetic semiconductors^ represent an important class of materials and struc- 
tures where ferromagnetism can be tuned by voltage^. This ability comes from the carrier- 
mediated character of the ferromagnetic interaction^^. The ferromagnetic ordered state in 
these systems appears due to mobile carriers interacting with stationary spins of magnetic 
impurities. To date, Curie temperatures as high as 40 K have been observed in a techno- 
logically important class of the Mn-doped III-V semiconductor structures'^. When the 
magnetic semiconductors become combined with the conventional field-effect layered struc- 
tures, the number of mobile carriers and the ferromagnetic interaction become tunable by the 
voltage 9 . This ability to externally control the properties of magnetic crystals with means 
other than the external magnetic field may have important device applications. A further 
step from magnetic semiconductor layers would be zero dimensional systems, quantum dots 
(QDs). Magnetic quantum dots can be viewed as nano-scale memory elements where infor- 
mation is stored in the form of magnetic polarization. Such a system may have important 
advantages compared to the conventional metal spin- valve memories: (1) small sizes and 
relatively small number of carriers and (2) voltage control of the number of electrons which 
was already demonstrated in many experiments for non-magnetic QDs 10,11 . Therefore it is 
important to develop the understanding of magnetic QDs with interacting carriers. 

Here we develop a theory of magnetic QDs with many carriers where Coulomb, ferro- 
magnetic, and single-particle energies contribute to the formation of the equilibrium state. 
Using the quasi-classical description, we show that a QD may be split into three phases 
with different physical properties. The geometrical sizes of these phases are determined by 
the Coulomb, ferromagnetic, and single-particle contributions to the chemical potential of 
a QD. For calculations, we employ the mean-field theory and the Boltzmann kinetic equa- 
tion. This approach becomes reduced to the Thomas- Fermi model at low temperatures. We 
should note that our approach ignores the discrete structure of single-particle spectrum of 
QD and is valid when electrons occupy at least several quantum levels. At the same time, 
this approach has an important advantage: it allows us to describe the Coulomb effects in 
relatively large QDs where, as it is shown below, the Coulomb interaction is very strong and 
significantly exceeds the ferromagnetic interaction and the kinetic energy of carriers. The 
hole-mediated ferromagnetism in quasi-two-dimensional (2D) systems is strongly anisotropic 
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due to the heavy hole-light hole splitting in the valence band. Therefore, the magnetic po- 
larization occurs predominantly in the growth direction. Then, two magnetic states of a QD 
with spin polarizations "up" and "down" can represent a single bit. 

Presently QDs and other nano-structures doped with magnetic (Mn) impurities attract a 
lot of attention. Among other studies, several recent theoretical papers investigate QDs and 
their electron and excitonic states in the presence of a single Mn io n 12 i 13 i 14 ' 15 . In particular, it 
was suggested in ref.— that a QD with a single Mn ion can act as a multi-qubit which can be 
controlled optically. Another direction of research describes the magnetic states and polarons 
in QDs with many Mn ions and one or several carrier o 16 ! 17 ! 18 ' 19 ' 20 ' 21 . Among the above 
publications, the paper— demonstrates that the Coulomb-interaction effects in few-electron 
QDs can determine a collective magnetic state of holes and Mn spins. Ferromagnetism and 
spin separation in digital layered structures and quantum wells were also studied in several 
experimental^* 2 ^ and theoretical publications 5,23,24 . 



II. MODEL 



As a model system, we are going to use a QD "made out of a 2D quantum well. Such a 
zero-dimensional system can be fabricated by etching and lithographical methods. Within 
the lithographical method, a QD can be defined, for example, by using the top metal gates 
(fig. QJi). The number of carriers in such a QD is a voltage tunable parameter. 

To describe a state of many carriers confined in a QD, we start from the local properties 
of the coupled hole-Mn system in a 2D quantum well. In our system, a mobile hole and 
Mn spins experience an exchange interaction: U exc = — (3/3 Y^i(Sz,dz)8{ r h — where 
(3 is the exchange interaction constant, and S z and j z are the z-components of the spin 
operators for a Mn spin and hole, respectively; and R, are the coordinates of hole and 
i-impurity, respectively. The above operator describes the interaction between Mn-spins and 
heavy holes and assumes a sufficiently large energy separation between the heavy-hole and 
light-hole quantized states in the valence band. The corresponding effective spin-dependent 
potential of a single hole has a form: 



U S pi n {j z ) 



3z ttQ 
| spim 



(1) 



3 




N 2D *10 12 (cm" 2 ) 

FIG. 1: (a) Model of a Mn-doped lithographically defined QD based on a semiconductor quantum 
well. The number of particles is controlled be the voltage applied between the top and back gates. 
The QD confinement for holes is given by a voltage difference AV = V\ — V 2 < 0. (b) Calculated 2D 
spin density as function of the total 2D density of holes at various temperatures; m^h = 0.38 tjiq, 
L = 70 A, iVo = 23 nm~ 3 , (3Nq = —2.2 eV, x e ff = 0.01. Curie temperature for the above 
parameters is about 6.4 K. 

where 

U° pin (N spin ) = -^x eff N J dz^(z)S Mn B SMn [^N spm ^ 2 (z)]. (2) 

Here N spin = N +3 / 2 — iV_ 3 / 2 is the net spin 2D density, Nj x are the 2D densities of the 
spin components, B$ Mn is the Brillouin function, Nq is the number of cations per unit cell, 
and Smu = 5/2. N 2D = N +3 / 2 + N_ 3 / 2 is the total 2D density in a quantum well. For the 
ground-state wave function in a quantum well, we will use ipa{z) = a/2/L sin7rz/L, where L 
is the quantum well width. The chemical potential for a 2D gas depends on U® in and N 2 d' 



U° ■ / U° ■ 2E 

H2d(T, N 2D , U spin ) = k B Tln[-cosh(-^) + \ cosh{-^) 2 + exp(—t) - 1], (3) 

Kb-L y kb± KbJ- 

where Ef(N 2 u) = nN 2 r)h 2 /rrihh- Now we calculate the spin polarization: 



^2D~ U s P in( N syin) 

N . = k B TD w 1 + e ^ 
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where D 2 d = mhh/i'xfo 2 )- The Zener ferromagnetic phase transition occurs when eq. 0] 
has a nonzero solution. Fig. shows the data for spin density N sp i n for a GaAs/AlGaAs 
quantum well with the following parameters: ra^h — 0.38 m , L = 70 A, N = 23 nm~ 3 , 
/3N = —2.2 eV, x e ff = 0.01. The above exchange parameter (3 is comparable to that 
used in other publications on magnetic semiconductors (see e.g.—). Since the exchange 
interaction is antiferromagnetic (/3 < 0), the case N spin > corresponds to the negative 
average polarization of Mn ions, Bs Mn < 0. Curie temperature can be analytically calculated 
in the high-density limit: k B T c = S(S + l)f3 2 x eff N m hh /(8nh 2 ) (k B T c < N tot /D 2D ). 

The in-plane potential in a lateral QD near its center can be approximated by the 
parabolic function: 



e4(r) = U + — - — r , (5) 

where e > is the electron charge and a; is a characteristic frequency of a confining potential. 
The potential Uq determines a depth of a lateral potential well. In a QD defined by metal 
gates (fig. ^,), the parameters Uq and to are functions of the gate voltages. In equilibrium, 
the carrier distribution function, which satisfies the Boltzmann equation, has a form: 



f^P^^z) P 2 /2m hh +e^r) + U SJ>in U z ,r)-„ ' ( 6 ) 

e k B T +1 

where r = (x, y) is the lateral position vector and p is the in-plane momentum. The self- 
consistent scalar potential of a hole is composed of two terms: 



e0(r) = e0 o (r) + U Co ui(r), (7) 

where Ucouiij) is the electrostatic potential induced by a non-uniform spatial distribution 
of carriers, n 2 z)(r) = n +3 / 2 (r) + n_ 3 / 2 (r). In addition, the distribution function (jHJ) depends 
on the spin of hole through the exchange interaction which is a function of the local spin 
density, n spin (r) = n +3 / 2 —n_ 3 /2 (see eqs. HI2|) . At the same time, the function n spin (r) itself is 
determined by the total local density of holes, n2_o(r), and is given by the numerical solution 
of eq. |U (see the data in fig. [TJd). Therefore, it is convenient to regard U sp i n as a function of 
n 2D , i-e. U spin [j z ,r] = U spin [j z ,n spin {n 2D ))} = U spin [j z ,n 2D (r)]. By integrating the function 
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(JHJ) over momenta we come to two non-local non-linear equations for the densities w±3/2(r). 
Then, these equations can be solved for the chemical potential and rewritten in the form 
resembling the central equation of the Thomas-Fermi model: 



2irh^ rij 

/i = e0 o (r) + U Coul (r) + (j z / \j z \)U^ n [n tot (r)] + k B T\n[e^^ - 1], Jz = ±3/2, (8) 

where j z = ±3/2 and U Co ui(r) = e 2 J ^ £eff( y°^\ r >_ r \ , where e e// (r' - r)| is an effective 
non-local dielectric constant of a system with metal gates^. We should also note that 
Ucoui{ r ) was written as a 2D integral and this it valid if the lateral size of a QD is greater 
than the quantum-well width L. 

In the system with the top gates closely located to a quantum well, the 2D integral in 
the equation for Ucoui( r ) is reduced to a linear function of n 2 ir^ and is given by a local 
flat-capacitor formula: 



Ucml(T) = 4 ^"^M (9) 

where d\ is the distance between the QD plane and the top gate; the distance to the back 
metal contact is assumed to be larger, i.e. d 2 3> di, e s is a dielectric constant of the 
semiconductor (e a = 12.5). The approximation Q has been successfully used in the past for 
description of several experiments on optical and electronic properties of modulated lateral 
structure o 2 ^ 2 ? . By using the local approximation for the self-consistent potential 0, we 
reduce eqs. El to coupled non-linear local equations which should be solved numerically. The 
total number of holes in a QD is determined by the chemical potential \x and the lateral- 
potential depth Uq. Electrostatics of the structure under study (fig. Hk) is similar to that 
studied in refs . 26 ' 27 ' 28 and we can use here the results of the above publications. If the barrier 
between the QD and back contact permits efficient tunnelling, the chemical potential in the 
QD coincides with the potential of the back contact (i.e. fi = 0). Simultaneously, the front 
barrier (usually made of AlGaAs) blocks tunnelling between the QD and the top gate. Also, 
if di <S d 2 , the potential U becomes close to eV\. 
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III. MAGNETIC PHASES IN QUANTUM DOTS 



Figs. 12131 and 0] show numerical calculations for the local spin densities in a circular QD 
with d\ = 300 A, Hon = 2 meV, and fi = 0; for the QD depth, we take U = —0.08, —1, 
and —2 eV. A QD with the minimum free energy is circularly symmetric and can be 
split into different phases. The total number of holes in a QD is given by an integral 
Ntot,QD = / n 2D {r)d 2 r. The corresponding N totj Q D for the above values of U are estimated 
as 23,3550, and 14200. In QDs with relatively small N to t,QD (hglH), the system is split into 
ferromagnetic (F) and paramagnetic (P) phases. In fig. |21 the carriers with spin j z = +3/2 
are pushed away from the center of QD, the total spin of holes is negative, and the Mn- 
subsystem has a positive magnetization. This situation corresponds to the antiferromagnetic 
hole-Mn coupling (/3 < 0). With increasing the total number of carriers (fig. EJ), one can 
see the formation of another stripe within the ferromagnetic phase. This stripe is located 
between the center region of a QD and the paramagnetic phase and the holes in this stripe 
are almost fully spin-polarized. With further increasing N to t,QD (fig- EI) and for relatively 
low temperatures, the formation of the ferromagnetic stripe (F2) with filly-polarized holes 
becomes evident. Simultaneously, the paramagnetic stripe becomes very narrow. Such a 
magnetic stricture of a QD can be understood by looking at the data in fig. ^p. At low 
temperatures, the function N spin {N 2 £,) becomes very close to the linear function N 2 d i n an 
extended interval of N 2D . For example, at T = 4 K, N spin (N 2 D) ~ N 2D for 0.2 * 10 12 < 
N 2 £> < 3 * 10 12 . In the above interval of N 2 d at T = 4 K, the hole subsystem is almost 
completely spin polarized. 

It is interesting to estimate different types of energies contributing to the formation of 
stripes. It is easy to see that the Coulomb energy in eqs. |H1 dominates the ferromagnetic and 
single-particle (kinetic) energies. The Coulomb energy Ucoui — 4:ire 2 dN 2 u/e s ~ 90 meV for 
N 2 d = 2 * 10 11 crrT 2 while the spin energy |£/° pi J ~ 2 meV for N spin = 2 * 10 11 cm~ 2 at 
T = 4 K. The single-particle kinetic energy under the similar conditions Eki n ~ ti 2 d/D 2 o ~ 
1.2 meV. For QDs with more carriers and higher n 2 D, the above energies become increased 
but the condition Ucoui » \Us pin \ ~ E kin remains. 

We can also obtain analytic solutions of eqs. |H1 under certain conditions. If both spin 
subsystems of holes (j z = ±3/2) form a degenerate Fermi gas, the last term in eqs. |H1 
becomes proportional to the Fermi energy, 2nj z /D 2D . Then, we can sum up the equations 
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FIG. 2: Calculated hole density as a function of the distance from the center of QD for a relatively 
small number of holes. The QD is divided into ferromagnetic (F) and paramagnetic phases (P). 
The dashed line represents the total 2D density. Inset: sketch of the QD structure. 

for j z = ±3/2. The resulting equation does not contain the spin energy Ug pin and can be 
solved analytically: 



This formula can be applied, for example, to large QDs in the spatial region of the fer- 
romagnetic phase Fl, < r < Rpi, where Rpi is the radius of the Fl phase (see 
fig. HJ). For this phase, spin densities can also be found analytically, by using the condi- 
tion n spin = n_ 3 /2 — n +3 / 2 = —N spint 

saturation^ 

where N spin>saturation is a positive constant 
equal to N spin at high N 2D in fig. Wp, for T = 4 K, N spintSaturation = 3.26 10 12 cm" 
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fig- lit)- The formula fllOD also describes the density distribution in the paramagnetic phase 
in the regions where the hole gas is degenerate (Ef = U2d/D2d > fcgT). For many other 
regimes, the spin densities should be found numerically. Since the Coulomb energy domi- 
nates the magnetic and kinetic terms, the total radius of a QD can be well estimated from 
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Distance, r (nm) 

FIG. 3: (Calculated hole density as a function of the distance from the center of QD for a larger 
number of holes. The QD is divided into ferromagnetic (Fl and F2) and paramagnetic regions (P). 
Inset: schematics of the QD structure. 

eq. [TUlby putting tl2d{Rqd) = 0. The resulting estimate Rqd ~ ^^\Uo\/[m hh u; 2 } is valid at 
low temperatures. 

Experimentally, the stripe structure of a QD can be observed, for example, by spatially- 
resolved optical spectroscopy^. In optical spectroscopy, a spatial resolution can be as small 
as 0.1 [im^. Simultaneously optical emission is sensitive to the spin-polarization of carriers. 
In an optical experiment, a ferromagnetic QD would be excited with weak nonpolarized 
illumination; the resulting local photoluminescence will be circularly polarized and reveal 
the formation of stripes with different magnetic structures. 

Optical methods can also be used to write a magnetic state of QD (bit: "up"-"down"). 
This may be done with circularly-polarized light. Polarized optical pumping can bring a QD 
into a required collective magnetic state with spins "up" or "down". In order to prepare a 
quantum dot in a required magnetic state, one can also use a magnetic field induced by an 
electric current driven through a metallic wire on the surface of a sample. 

Another method to observe the magnetic phases in nano-structures is the electrical- 
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Distance, r (nm) 

FIG. 4: Calculated hole density as a function of the distance from the center of QD for a large 
number of holes at low temperature. The QD is divided into ferromagnetic (Fl and F2) and 
paramagnetic phases (P). Inset: schematics of the QD structure. 

capacitance spectroscopy^*^ which was successfully applied to observe, for example, com- 
pressible and incompressible stripes in electron quantum wires in the regime of the quantum 
Hall effec t 27 i 28 . The capacitance spectroscopy has been also applied to lateral and self- 
assembled quantum dots (see e.g.—). For the nano-structures with relatively large sizes 
considered in this paper, the signature of the ferromagnetic phase in the capacitance spec- 
tra is expected to be relatively weak because of the inequality Ucoul » \Ug pin \. However, 
the formation of the ferromagnetic phase can be recognized from a critical behavior of the 
capacitance spectrum as a function of temperature and voltage. 

IV. CAPACITANCE SPECTROSCOPY AND MAGNETIC PHASES IN QUAN- 
TUM WIRES 

As an example, we consider here quantum wires in a structure with the interdigitated 
metal gate (see inset in fig. |5J). In such a system, alternating voltages are applied to the 
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metal strips located on the surface of a quantum well. Similar structures were studied 
experimentally in ref.— . We can calculate the capacitance of a wire as a derivative C(V\) = 
e(dN totj Q W I dVi)l ~ e 2 (dN totj Q W /dU )l, where N ta t,QW is the linear density of carriers in 
a quantum wire and / is the length of a wire in the in-plane y-direction. In the local 
approximation for the Coulomb potential (eq. |HJ), the problems of quantum dot and wire 
become similar. Fig. |3] shows the capacitance of a quantum wire as a function of voltage 
for two temperatures, just below and above the Curie temperature Tq = 6.4 K. One can 
see in fig. El that at low temperatures (T < T c ) the capacitance becomes increased starting 
from a critical voltage Uo/e ~ 0.02 V. This voltage corresponds to the minimum 2D density 
(~ 5 * 10 10 cm -2 at T = 4 K; see fig. \^p) which is necessary to obtain the ferromagnetic 
phase stripe in the center of nanowire. Starting from this voltage, the central region of a 
nanowire contains a ferromagnetic stripe. The capacitance of a partially ferromagnetic wire 
becomes increased since the spin interaction makes a lateral potential well a little deeper 
and a wire can accommodate more carriers at a given voltage. If temperature increases just 
by 3 K, the ferromagnetic stripe vanishes and capacitance becomes reduced. This peculiar 
temperature behavior for Uq > 0.02 eV can be taken as an evidence for the ferromagnetic 
phase. 



V. CONCLUSIONS 



In this paper, we studied quantum dots and wires with many interacting carriers within 
the quasi-classical approach. The strongest interaction in quantum dots with a relatively 
weak confinement and a large number of carriers comes from the Coulomb forces. However, 
a weaker ferromagnetic interaction determines the spin structure of a large quantum dot. 
Depending on the parameters, a quantum dot can be split into three phases. 
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FIG. 5: Calculated capacitance of a nanowire in the system with interdigitated metal gate (see 
inset on the left hand side) for two temperatures, 4 and 7 K. The inset on the right hand side 
shows a top view of a nanowire with two phases. 
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